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Optical flares from the faint mid-dM star 2MASS J00453912+4140395 
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We present B and V light curves of a large stellar flare obtained with the Wide Field Camera at the Isaac Newton 2.5-m 
telescope (La Palma). The source object is a faint {my = 21.38) foreground star in the field of the Andromeda galaxy, 
with its most probable spectral type being dM4. We provide an estimate of the total flare energy in the optical range and 
find it to be of the order of 10 35 erg. The cooling phase of the large flare shows three additional weak flare-like events, 
which we interpret as results of a triggering mechanism also observed on the Sun during large coronal mass ejections. 
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1 Introduction 

Flares are known as sudden and violent events releasing 
magnetic energy and hot plasma from the stellar atmospheres. 
We observe them on magnetically active stars and, much 
more closely, on the Sun. Electromagnetic radiation is emit- 
ted across the entire spectrum, from radio waves through the 
optical range to X-rays and 7 rays. The total energy released 
during a typical solar flare is of the order of 10 30 erg, while 
the largest solar two-ribbon flares can emit up to 10 32 erg. 

According to the accepted model of solar flares, in the 
upper atmosphere, between oppositely oriented magnetic 
field lines, a current sheet forms and magnetic reconnection 
takes place, which results in acceleration of particles, and 
thus produces electromagnetic radiation and plasma heat- 
ing. On the basis of the solar paradigm one could simply 
expect that modelling stellar flares is just a question of scal- 
ing (e.g., as a function of the released energy, size, dura- 
tion, etc.). However, besides general similarities, stellar flare 
observations also unveil problems that cannot be explained 
with the extended canonical solar flare model, arising from 
the different spectral distribution of the emitted energy, the 
role of age and spectral type of the host star, the multiform 
magnetic field topologies on stars, the tidal forces in active 
binaries, etc. The only way we can come closer to draw- 
ing up similarities and differences between solar and stellar 
flares is to continuously broaden the ensemble of analyzed 
flares (solar and stellar as well), thus making an observa- 
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tional base of different types of flaring activities available 
for theoretical purposes. This is how solar flare observations 
can help in the correct interpretation of stellar flares, and re- 
versely. 

Flare activity in cool stars is a very common phenomenon. 
Flares in stars with spectral types earlier than M are ob- 
served mainly in UV and X-rays and optical flares (like 
the most energetic white light flares in the Sun; Hudson, 
Wolfson, & Metcalf 2006) are rare. However, in the less lu- 
minous low-mass dM stars optical flares often occur. It is 
also known that red dwarfs from » 0.3M Q down to the hy- 
drogen burning limit can release a considerable amount of 
energy via flare eruptions, even when, in most cases, these 
objects show very low magnetic activity in their quiescent 
state. In this paper we present photometric observations of a 
large optical flare event observed on 2000 September 25 that 
is associated with one such low-mass dM stars. In Sect.|2]a 
short description of the data is given, and in Sect. [3] we sum- 
marize the stellar properties of the host star. In Sect. [4] we 
give an estimate of the energy released by the large flare 
and, finally, in Sect [5] an additional two weaker short-term 
flares are identified and a plausible scenario is set up for the 
three weak post-flare events that followed the large flare. 

2 Data 

Our time-series data were collected between 1999-2003 within 
the course of a variability survey in in the North-Eastern 
quadrant of the Andromeda Galaxy (M31). The observa- 
tions were acquired with the Wide Field Camera at the Isaac 
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Table 1 Catalog data on the observed target 





INT WFC 


DIRECT 


2MASS 


object ID 


M31J00453912+4140395 


D31J004539. 1+414039.5 


J0045391 1+4140396 


a (J2000) 


0045 39.12 


0045 39.11 


0045 39.12 


5 (J2000) 


+4140 39.5 


+414039.46 


+4140 39.7 


B 


23.18 






OB 


0.11 






V 


21.38 


21.62 




a v 


0.03 


0.07 




I 




18.18 




Ol 




0.07 




J 






16.35 


crj 






0.11 


H 






15.98 


on 






0.16 


K 






15.63 


OK 






0.18 




30 w F— — J* 

LL-. I , C 1 

Fig. 1 Finding chart from the M3 1 field with our target 
(INT WFC image). 

Newton 2.5-m telescope (INT) in La Palma (Spain). The 
program is described in detail in Vilardell, Ribas, & Jordi 
(2006). In brief, the observations, taken during the course 
of 21 nights, were reduced and then analyzed by means of 
the so-called Difference Image Analysis (DIA) technique, 
which is especially tailored to detect variable objects. The 
resulting data are high precision (~ 0.01 mag) photometric 
light curves in the Johnson B and V passbands. The object 
studied here, whichhas identifier M31J00453912+4140395 
in the Vilardell et al. (2006) catalog, was flagged during the 
course of the analysis as an object with a large brighten- 
ing following the characteristic light curve shape of a stellar 
flare. The finding chart of the object, also cross-identified 
as 2MASS J00453912+4140395, is shown in Fig.Q] which 
contains only a small portion of the full 32'x32'field of 



view of the WFC. The B and V observations resulting from 
the DIA analysis and subsequent calibration to the standard 
system are plotted in Fig.|2l where the large flare event (A) 
is zoomed alongside with two other weaker flares (B, C). 

3 Stellar properties 

TableHJsummarizes the quiescent B and V magnitudes from 
our INT data, as well as V and I magnitudes from the DI- 
RECT catalogue (Mochejska et al. 2001), and J, H and 
K magnitudes from the 2MASS catalogue (Skrutskie et al. 
2006). A self-consistent analysis using the observed object 
color indices and the extinction law of Drimmel, Cabrera- 
Lavers, & Lopez-Corredoira (2003) indicates a color excess 
of E(B — V) = 0.05 mag, and, adopting a ratio of 3.1 be- 
tween Ay and E(B — V), leads to an absolute extinction of 
0.16 mag. The unreddened color indices are most compat- 
ible with a spectral type M4 dwarf with My = 13.8 mag, 
log L/L Q « -2.5, M « 0.16 M , and R w 0.2 R ac- 
cording to the tabulations from Baraffe & Chabrier (1996), 
Baraffe et al. (1998), and Bessell, Castelli, & Plez (1998). 
The distance resulting from the analysis is of roughly 300 pc. 

KPNO plates taken from the M31 field (NO AO Science 
Archive, Aladin Sky Atlas, http://aladin.u-strasbg.fr/) also 
agree with a mid-M or later spectral type classification since 
the object is best seen in the /-band KPNO plates, but also 
observable in Ha and R. However, in the U, B, and V band 
plates our target is indistinguishable from the background. 

A first step towards evaluating the level of magnetic ac- 
tivity of the studied target is to investigate possible cyclic 
(rotational) variability induced by surface inhomogeneities 
(e.g., cool spots). No convincing sign of such modulation 
was found from our Fourier analysis. Three possible expla- 
nations can be put forward. First, spot coverage could be 
too small to observe any rotational modulation (this is the 
case when observing our Sun as a star in optical photomet- 
ric bands). Alternatively, the star could be heavily spotted 
and yet show small or no modulation if spots are evenly dis- 
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Fig. 2 Light curves in Johnson B and V bands of 2MASS J00453912+4140395 taken with the WFC of INT in La Palma, 
between 1999-2003. Three individual flare events (flares A, B and C) are marked and zoomed. 



tributed. If we compare the mean error of the measurements 
in B is 0.08 mag with a scatter of 0.11 mag, i.e., there is 
still room for some low-level modulation, although the V 
band does not seem to show such behaviour with a scatter 
of 0.03 mag and a 0.04 mag formal error. The third possible 
scenario is a low inclination (i « 0°) of the rotation axis 
making it not possible to observe any rotational variability. 
However, in this latter case one could expect variability on 
a longer timescale of a few years, because of changes in 
the overall spot coverage according to solar-like spot activ- 
ity cycle (see e.g. Olah, Kollath, & Strassmeier 2000). No 



such long-term modulation is seen during the 5-year long 
observing season and 2MASS J00453912+4140395 seems 
to be like the vast majority of cool M dwarfs with rotational 
variability less than 1 — 2% in the visible. 

4 Flare energy estimation 

Flare statistics of red dwarf stars show that stellar flares can 
generally be divided into two subgroups: a group consist- 
ing of relatively small, short (w 10 3 s) impulsive flares, and 
a more energetic group releasing at least 10 32 erg and last- 
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ing w 10 4 s. The latter are often related to solar two-ribbon 
flares which are associated to filament eruptions and coronal 
mass ejections (CMEs). The basic properties of flare A ob- 
served on 2000 September 25 (beginning at HJD 245 1 8 1 3 .45), 
such as the long duration reaching Ats ~ 3 x 10 4 s and the 
large amplitudes of AB =2.69 mag and AV =1.23 mag 
clearly class this event among the more energetic group. 

Flare light curves usually consist of a rapid rise followed 
by a slower, monotonic decay. However, our poorly covered 
rising phase permits only a rough estimate of the physical 
properties. For the estimation of the flare energy first we 
derived intensity from the magnitude values: 

!o+f 



= 10" 



(1) 



where 7o+/ and 1$ are the intensity values of the flaring 
and the quiescent stellar surfaces, respectively, in one of the 
observed bands. The relative flare energy is then obtained 
by integrating the flare intensity over the flare duration: 



£ 



(2) 



The quiescent stellar fluxes in different bands are es- 
timated assuming a simple black body energy distribution 
for a dM4 star with T cS = 3100 K and R « 0.2 R (see 
Sect. [3]) from the equation 



4TrR 2 T(\)S B ,v(\)d\, 



(3) 



where .F(A) is the power function and Sb,v(A) is the trans- 
mission function for a given passband. Finally, the total inte- 
grated flare energy is calculated by multiplying the relative 
flare energy by the quiescent stellar flux: 

E f = SfF*. (4) 

Since our flare light curve is poorly covered, the abrupt 
rising phase with its real peak value may be estimated as 
even 1.5-2 mag brighter in B and ral.3 mag brighter in V 
than the measured maxima. This assumption corresponds to 
a photometric flare temperature of « 1.5 x 10 4 K at peak 
(cf., e.g., Ishida 1990; Ishida et al. 1991, and references 
therein). According to this, the observed and dereddened 
B — V peak of 0.29 mag is just an upper limit of a more 
reliable B — V « —0.16 mag. Calculated luminosities for 
the quiescent star, and the total flare energy for both the min- 
imum fit and the more realistic assumption are summarized 
in Table|2] together with the ratio of the flare energy and the 
quiescent stellar luminosity, often called equivalent flare du- 
ration, i.e., the time interval in which the star would radiate 
as much energy as the flare itself. The ratio between the flare 
energies Eb/E\/ — 2.01 (or 1.83 for the more realistic es- 
timation) is of the order of the statistical value of 1.60^032 
from Lacy, Moffett, & Evans (1976). Using the empirical 
correlation 

E v = 1.2 ± 0.08E B (5) 

we obtain an estimation of 4.27 (7.61) 10 34 erg for Ejj, 
From this we can estimate the total flare energy released in 



Table 2 Quiescent stellar flux and flare energy estima- 
tions. A minimum flare energy is estimated both from a min- 
imum fit of the observed data and also by a (more realistic) 
reconstruction of the rising phase (values in parenthesis). 



band 


quiescent flux 


flare energy 


equivalent duration 




10 29 [erg/s] 


10 34 [erg] 


lh] 


B 


5.09 


3.56 (6.34) 


19.4 (34.6) 


V 


12.76 


1.77 (3.47) 


3.9 (7.6) 



the optical range as being of the order of a few times 10 35 
erg, and resulting in a value of ~ 10 36 erg for the bolometric 
flare energy (cf. the review of Pettersen 1989). This rough 
but still realistic estimation (cf. the example in Pagano et al. 
1997 with ~ 3 times more brightness in U) indicates that the 
2000 September 25 flare is among the most energetic stellar 
flares ever observed. 

5 Weak flaring and post-flare events 

There are statistical evidences that flare-like transient phe- 
nomena, often called microfiares, are almost continuously 
present on M-type dwarfs. These events are thought to be 
one of the major sources of chromospheric and coronal heat- 
ing. We see them (sometimes as nanofiares) also on the 
Sun (e.g., Lin et al. 1984), however, because of being rela- 
tively weak and short-term events most of them remain un- 
resolved. 

For filtering out such short-duration events from the back- 
ground of our photometric data we searched for small am- 
plitude peaks occurred simultaneously in both photometric 
colours. We applied different filtering algorithms but finally 
a visual inspection proved to be most reliable and efficient. 
Two such events are apparent. One is just one day after flare 
A at HJD 2451814.55 (marked with B in Fig.© and another 
one at HJD 2452203.45 (flare C in Fig. Hi. There are some 
additional albeit less convincing ones, which are only 2-cr 
away from the photometric background, and therefore were 
disregarded. Following the method described in Sect. |4j we 
estimate the total energy of flares B and C to be from a few 
times 10 32 erg up to 10 33 erg in B, which are characteristic 
of short-term impulsive flares rather than microfiares. 

When analyzing the decay phase of the light curve of 
flare A, three other short term increases can be identified. 
Since those events occurred simultaneously in both colours, 
we assume them to be an additional three weak flares: one at 
HJD 245 1 8 1 3 .54 (D in FigO, another one at HJD 245 1 8 1 3 .67 
(E) and a third one just at the end of the large flare light 
curve at HJD 245 1 8 1 3 .74 (F), all of them lasting about At w 
8 x 10 3 s (in the case of flare F, the simultaneous rise in B 
and V together is considered as just the beginning of the 
flare event and its duration is extrapolated). Again, we esti- 
mate the total flare energies as 2 — 6 x 10 31 erg in B and 
about a half of this value in V. The occurrence of such weak 
post-flare events in the descending phase of a large long du- 
ration flare is a reminder of the mechanism described first by 
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Fig. 3 The large flare of 2000 September 25 (flare A in Fig. |2]l with three post-flare events (D,E and F). Left panel shows 
the B (dots) and V (circles) light curves together. In the right panel the zoomed V curve is plotted alone. 



Attrill et al. 2007, which probably occurs on the Sun during 
large CMEs. After energetic solar flares are associated with 
CMEs, the moving footpoints of the blowing up magnetic 
loop crosstalk with opposite polarity flux ropes from the 
randomly distributed magnetic carpet, or other favourably 
oriented magnetic concentrations (active regions), thus form- 
ing current sheets and triggering other (micro)flare events 
from time to time. This seems a plausible interpretation for 
the three weak flare-like events (D, E and F) observed dur- 
ing the decay phase of flare A. 
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